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Abstract 

The helicity modulus for a fluctuating type II superconductor is computed 
within the elastic medium approximation, as a probe of superconducting phase 
coherence and the Meissner effect in the mixed state. We argue that at the 
vortex line lattice melting transition, there remains superconducting coher- 
ence parallel to the applied magnetic field, provided the vortex line liquid 

retains a finite shear modulus at finite wavevector. 
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In the high temperature superconductors, fluctuation effects are believed to be impor- 
tant over a wide region of the H — T phase diagram [|1]-|IJ. Recently there has been much 
controversy concerning the effect of vortex line fluctuations on long range order of the su- 
perconducting wavefunction in the mixed state 0-0]. Part of this controversy has concerned 
the definition of the proper gauge invariant correlation function. In this work we reconsider 
the question of superconducting coherence by considering instead the behavior of the helic- 
ity modulus || of the superconductor, which has proven a valuable criteria for coherence 
in related superfluid []10| and Josephson array [111 models. As will be seen below, it is an 



intrinsically gauge invariant quantity. Furthermore, the helicity modulus is equivalent to the 
linear response coefficient between an applied perturbation in magnetic field and the result- 
ing supercurrent. Hence it more directly probes one of the most characteristic properties of 
a superconductor, the partial Meissner effect of the mixed state. 

Working in the elastic medium approximation, we compute the helicity modulus to low- 
est order in the fluctuation of vortex lines. We show that vortex line fluctuations have a 
dramatically different effect on the helicity modulus parallel versus transverse to the ap- 
plied magnetic field. For the parallel case, the helicity modulus gives total screening as 
in the Meissner effect. We argue that this total screening is unaffected by the vortex line 
lattice melting transition, provided the vortex liquid retains a finite shear modulus at finite 
wavevector. In the limit of an extreme type II superconductor with A — > oo, our results ex- 
plain recent numerical simulations by Li and Teitel |L^] which show clearly the persistence 



of phase coherence parallel to the applied field, well into the vortex liquid state. Similar 
results were obtained by Feigel'man et al. fI3 |, working with a related 2d boson model. We 



show a simple relation between their results and ours, which can be expressed in terms of 
the "winding number" of vortex lines. For simplicity, we carry out our calculation for an 
isotropic superconductor; the extension to the uniaxial anisotropic case is straightforward. 
We work within the London approximation, which should be valid provided one is not too 
close to H c2 . 



The Landau-Ginzburg Helmholtz free energy [Jlj] for an isotropic uniform superconduc 



tor, within the London approximation of constant wavefunction amplitude, can be written 

as, 

H = ^J J d 3 r{\V9 - A| 2 + A 2 |V x A| 2 } (1) 

where 9 is the phase of the superconducting wavefunction, A is the "bare" magnetic penetra- 
tion length, Jq = 0q/167t 3 A 2 with </> the flux quantum, and (0o/2vr)A is the magnetic vector 
potential. V x A = 27rf where f(r) = B(r)/0 o is the local density of magnetic flux quanta. 
The partition function Z is computed averaging over independently fluctuating 9 and 
A, subject to the constraint that (f(r)) = fz for a uniform average magnetic induction Bz. 
In evaluating Eq.([|), the integration is to be cut off at the core of a vortex line, so that the 
energy stays finite. 

In terms of the supervelocity v = V# — A, and its Fourier transform y q = J d 3 r e iq ' r v(r), 
the helicity modulus is defined [|I(J as the linear response coefficient between induced super- 
current Jqv and an applied twist in phase. If we take v 9 — > v q + 5v q fx in Eq.(ff]), then the 
helicity modulus in direction p, is, 



T,(q) 



d 2 T 



= J o fl - 7?^(V-<zm) \ ( 2 ) 
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where T = —T In Z is the total free energy, and V = L Z L\ is the system volume. Because of 
the symmetry with which V9 and A enter v, equivalently gives the induced supercurrent 
that flows in response to an applied perturbation in magnetic field, given by the vector 
potential SA q = 8v q pL. In evaluating T M (q), the physically relevant case is the limit — > 0. 
This follows from the convention of Baym |16l , where to describe a system with a current 
flowing in the direction jl, the appropriate thermodynamic limit is to take the system size 
L M — > oo first, followed by L u — > oo for the directions v _L jl (equivalently, q ■ 5A q = in 
the London gauge). 

Defining the vortex line density n by V x W = 27rn, one can write an arbitrary config- 
uration v q in gauge invariant form as, 



X is a smooth function which gives the longitudinal part of v q , while the transverse part of 
Vg is determined by V x v = 27r(n — f). Substituting Eq.(^) into the Hamiltonian (|IJ), and 
decoupling the n q and f q degrees of freedom by completing the square in f g , we get 

2lT 2 J ^ 2 1 , 2 2\ rr src ^ 



n 



E [?X q X- q + ? (1 + A 2 g 2 )5f 9 • 5i_ q + ITW n ? ■ n_ 9 | (4) 



where 8f q = f 5 — f? is the fluctuation of the magnetic flux density away from the value 
fg = n g /(l + X 2 q 2 ), which minimizes the Hamiltonian for a given vortex configuration n q . 
The partition sum is now an average over all smooth functions Xi ^f<? such that q-^f^ = 
(so that V • B = 0), and all vortex configurations n. The vortex line interaction in Eq.(|4]) is 
just the familiar London result [|14],[17] • 



The helicity modulus is evaluated by substituting Eq.(Q) into Eq.(|2|), and using the 
Hamiltonian (|j) to evaluate the averages over \q an d 5f q . Taking the limit — > one gets, 



I- T / x _ J o A V J\ _ 47T 2 J A 2 (ft x q)q(A x q)^(w ga w_^) 

mn „ W - nm i + ^ < ^ x + A v 



(5) 



For a superfluid or spin model P-pll phase coherence is indicated by a non-vanishing 
Y M in the limit q — > 0. For the superconductor however, the gauge field A is free to 
adjust itself to screen out the applied phase twist (or perturbation 5 A), and so even in 
the superconducting state T^(q — > 0) ~ q 2 , as seen in Eq.([|) above Hl8f . In fact, if no 
vortex lines are present (n 9 = 0), Eq.([5|) just gives the familiar total screening response 
of the Meissner state |Tj|. With the presence of vortex lines in the mixed state, we can 



generalize the form of the Meissner response, by defining a renormalized coupling (JX 2 )r 
and penetration length Xr such that 



_ (JX 2 ) R q 2 
1 + X 2 R q 



v ~ ^ nq , a s 

llm , T M = 1 I ^2 ( 6 ) 



where Xr and (JX)r may depend on the direction q. Thus to examine superconductivity it 
is necessary to consider the form of T M at small but finite q. 

At high T, one can make a hydrodynamic approximation |H| and average over n(r) as if 



it was a continuous function, subject to the constraint that vorticity is conserved q • n q = 



Using the Hamiltonian (|J) one gets T /U (q) = as expected. At low T, one can evaluate 
the vortex line correlations using the elastic approximation for small fluctuations about a 
vortex line lattice. If Ui(z) is the transverse displacement of vortex line i at height z from 
its position R, in the vortex line lattice, then n(r_i_, z) = J2i <5( r _i_ — Ri — u i{ z ))( z + dvn/dz). 
To evaluate Eq. (^) to lowest order in T, it is only necessary to consider the expansion of n q 
to linear order in Uj. For small q > we have, 

n y = if(q- u q z - q z u q ) (7) 

where rj_, Rj, and Uj lie in the xy plane, q z and q_i_ are the components of q parallel 
and perpendicular to z, and u q = fY.il dz e^ qzZ+Ci± ' R ^u i (z) . Correlations of u q may be 
evaluated using the elastic Hamiltonian, as derived by Brandt ||17|| , 



Hei = Tfrp H {( c 44<?z + c u q 2 A_)u qL u^ qL + [c M ql + c m q 2 ± )u qT u^ qT } (8) 

where u q L and u q T are the components of u q parallel and transverse to q_i_, and C44(q), Cn(q), 
and C66(q) are the tilt, compression, and shear moduli respectively. 

Substituting Eq.(|7|) into Eq.(||), and evaluating the displacement correlations using Ti eh 
we find for perpendicular and parallel responses, 

in T ^ = £55, TTAV I " 4.(l + AV)(c 44 ^ Cll gi)) (9) 

Km T.(q) = lim {l - 47r (i + AV) + ^ ± ) } ' 

For the transverse response T x (q) there are two cases to consider: (z) q = qz, and (n) 
q = qy. In (z) the perturbation 5A q gives a magnetic induction along y, oscillating in the 
z direction. It is thus a tilt modulation of the original induction Bz. In (ii), the perturba- 
tion gives a magnetic induction along z, which oscillates along y; it is thus a compression 
modulation of Bz. Accordingly, Eq.(|9]) shows that in (i) T x depends on C44, while in (ii) T x 
depends on c\\. We consider in detail case (i). A comparison of Eq.(||) with Eq.(||) shows 
that (J\ 2 )r is determined by C44(q = 0), while Ar is determined by dc^O) / dq 2 . Using the 



result of Brandt |T7j 



B 2 

C44 = — 

4:71 
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1 + A 2 g 2 \dB A 



(11) 



we find 



J A 2 difj.' A 2 diii" 1 J 

For an isotropic system, the factor dH±/dB± in C44 above, where the derivative is evaluated 

at the average magnetic induction Bz, is equal to the more familiar Hj B. The renormaliza- 

tion factor for the coupling (JX 2 )r has a simple physical interpretation. Since the induced 

magnetic induction is determined from Maxwell's equations as = — (T x (q)/JoA 2 g 2 )5A, 

Eq.(|l2D results in a fraction dB±/dH± of the perturbation 5Hy = V x 5 A penetrating the 

superconductor, while the remainder is screened out as in the Meissner effect. 

We now consider the parallel response T 2 of Eq. fllQ|) . As long as the shear modulus c 66 is 
finite in the limit q z — > 0, the term in Eq.([U]) due to vortex line fluctuations vanishes, and 
one has total screening of the perturbation as in the Meissner state. If cqq is identically zero, 
T 2 (q_|_) has exactly the same form as T x (qz), with the greatly reduced coupling l — dB±/H±. 
A related analysis of order parameter correlations || has led Ikeda et al. to conclude 
that vortex line lattice melting (where C66(0) — > 0) should be accompanied by a dramatic 
reduction in phase correlations along the direction of the magnetic field. However a careful 
analysis of Eq.(|i"0|) shows that T 2 continues to be unaffected by vortex line fluctuations 
provided that Cggfe, q_i_) does not vanish as fast (or faster than) g 2 , as q z — > for finite q_i_. 

One may question the application of elastic theory once the vortex line lattice has melted. 



A possible justification has been given by Marchetti and Nelson who show that a hexatic 
vortex line liquid may be described by an elastic theory in which one includes free dislocation 
loops. Averaging over dislocations, they find that the elastic moduli en and C44 are largely 
unchanged, however the shear modulus is renormalized to Cgefe = 0, qjj ~ q\. Thus, 
while c 66 (q = 0) = as expected for a liquid, c 66 remains finite for finite q^. Hence 
Eq. ( |I0D continues to give total screening of the perturbation upon melting of the vortex line 
lattice, ie. superconducting coherence persists parallel to the applied field for some finite 
temperature range into the vortex liquid state. 



Continuing the expansion as in Eq.(^) to next order in the displacements u q , we find a 



correction only to \r of order A|,/A 2 ~ (3.8T/7rJ )y B/(f) (using B ~ 0.2H c2 ). Evaluating 
at the vortex line lattice melting temperature, which we find to be T M ~ 1.7 'c^ttJq J q>o/ B 



(where cl ~ 0.15 is the Lindemann parameter), we find a small correction to X 2 R of order 
15%. Thus the conclusions above from the lowest order expansion continue to hold ||22j| . 

If one considers the above calculation in the limit of an extreme type II superconductor 
where A — > oo, and all fluctuations of the gauge field A are frozen out, we return to the 
case analogous to an ordinary superfluid. Taking the A — > oo limit in Eqs.(|9]— 11), we find 



that T z (q — > 0) = Jo is finite, and hence the system has phase coherence in the z direction, 
even in the vortex line liquid state (provided > for qj_ ^ 0). however vanishes at 
all temperatures, as B = H when A — > oo. This explains the recent numerical results of 
Li and Teitel 0] in a lattice A — > oo model. There T z was found to vanish at a T cz , well 



into the vortex line liquid state, while T x vanished at a much lower T c ±, where the vortex 
line lattice melted. The finite T x for this model at low T is due to the effects of pinning 
introduced by the discrete numerical mesh, which creates a finite energy barrier to small q 
elastic distortions. This effectively adds a q independent constant to the denominator of the 
second term in Eq.(^), so that as q — > 0, = Jo- Only when the vortex line lattice melts 
will thermal fluctuations dominate over pinning, and one recovers Eq.(|i2|) with the resulting 
T x = 0. 

Much work has been done using an analogy between fluctuating vortex lines, and the 
imaginary time world lines of two dimensional bosons |I|Jl3|. Feigel'man et al. ]13| have used 



this analogy to argue that in the large A limit, there will exist a boson normal fluid phase 
intermediate between the boson lattice and the boson superfluid phases. They predict that 
the corresponding phase of the superconductor is characterized by coherence parallel to the 
applied magnetic field, but not transverse to it. It is interesting to examine this analogy 
within the above elastic approximation. A convenient expression for the 2c? boson superfluid 
density has been given by Ceperley and Pollock E3] in terms of the "winding number" W 



of boson world lines, p s = mT(W 2 )/2h . If there are only magnetic field induced vortex 



lines fluctuating in a directed fashion [|4j] (ie. a single valued displacement Ui(z)), 

(W 2 ) = ^(| EWL.) - u,(0)]| 2 ) = -L« =0 ■ n^ =0 ) (13) 



where n^ =0 is the average vortex line density transverse to the average magnetic induction 
Bz. Translating [I]] from 2d bosons to vortex lines (h/Tb oson — > L z , h —>■ T vortex , tti — > e\ ~ 
7rJ the single vortex line tension), and evaluating (FT 2 ) within the elastic approximation, 



one finds [25 



p s = lim lim { J* T + J* I . (14) 

9±^o qz ^o 2 { c AA q 2 z + c m qj_ c AA q 2 z + c u q z ± 



The second term above always vanishes when one takes q z — > first, as en is always finite. 
The first term is just the same factor as appears in Eq.(pl]) for T 2 . Hence in the vortex line 
liquid, if c 66 (q_i_ ^ 0) > 0, we have both total Meisnner screening of perturbations 5A q± z, 
and p s = 0, consistent with the predictions of Feigel'man et al. 

When c 66 vanishes identically, Eq.(|l4]) gives p s = eif 2 /2c AA {$). This result also follows 
from a direct evaluation of Eq.(|l"3"D within a hydrodynamic approximation ||19|| . If the 2d 
boson normal to superfluid transition is of the Kosterlitz-Thouless type, then the universal 
jump in p s at the transition may be written |23J as (W?) = A/tc. This gives a transition 



temperature for the vortex lines, Tkt — {4>o/^ 2 L z )dH±/dB±. For H ~ H&, this result 



gives Tkt ~ in good agreement with an earlier estimate by Nelson | 26fl . At larger B, 
where dH±/dB± ~ 1, Tkt is independent of B. For sufficiently large L z however, this KT 
transition is presumably preempted by the transition to the hexatic vortex line liquid, in 
which c 66 (q_i_) > 0. 
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